Competing magnetic fluctuations in 7— iron due to peculiarities of its electronic 

structure 
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Applying the local density and dynamical mean field approximations to paramagnetic 7-iron, 
we revisit a problem of theoretical description of its magnetic properties in a wide temperature 
range. We show that contrary to a-iron, the frequency dependence of the electronic self-energy 
has a quasiparticle form for both, t^ g and e g states. The momentum dependence analysis of static 
magnetic susceptibility shows maximum near the wave vector qx = (27r, 0, 0), in agreement with the 
experimentally observed antiferromagnetic state in 7-iron precipitates. This state is however found 
to closely compete with the states, characterized by magnetic wave vectors along the directions 
qx-(7r,7r,7r)-(37r/2,37r/2,0). From the analysis of the uniform magnetic susceptibility we find, that 
contrary to a-iron, the Curie- Weiss law is not fulfilled in a broad temperature range, although the 
inverse susceptibility is nearly linear in the moderate-temperature region (1200-1500 K). This non- 
linearity of the inverse uniform magnetic susceptibility in broader temperature range is due to the 
density of states peak located in the vicinity of the Fermi level. Based on the obtained momentum 
dependence of the susceptibility in the paramagnetic phase the exchange integrals are estimated. 
The obtained results allow us to ascertain the size of the correlation effects in 7-iron as weak to 
moderate and give evidence of close relation of 7-iron to the other weak itinerant antiferromagnets. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

The problem of magnetism of iron attracts a lot of at- 
tention till now. Pure a-iron has body centered cubic 
crystal lattice and it is ferromagnetic at temperatures 
below Curie temperature 1043 K 1 3 . In the tempera- 
ture range between 1043 and 1183 K a- iron is param- 
agnetic. This most studied allotrope of iron becomes, 
however, unstable above 1183 K because of the struc- 
tural phase transition to the 7-phase 1,4 , which has a face 
centered cubic crystal structure 2 ' 3 . The theory of the 
a- j structural transition is still under development. Re- 
cent investigations 5 7 have shown important role of mag- 
netic correlations for this transition. These observations 
are supported by the results, indicating presence of local 
magnetic moments in a-iron even above the magnetic 
transition temperature 8 . 

In view of these observations, understanding of mag- 
netic properties of 7-iron seems to be necessary. Exper- 
imentally, the temperature dependence of inverse mag- 
netic susceptibility in 7 phase has very weak slope, which 
cannot be determined to a good accuracy because of large 
spread of experimental data (see Refs. 9,10 and refer- 
ences therein). Curie- Weiss temperature extracted from 
a fit to experimental data is negative, Ocw — —3451 K, 
and the corresponding magnetic moment is about /icw = 
7.47/iB 10 - Therefore, magnetic properties of 7-iron are 
very different from those of a-iron, for which the Curie- 
Weiss temperature Ocw — 1093 K and the magnetic mo- 
ment /icw = 3.13/iB 9 - The particles of the size between 
60 A and 700 A in precipitates of 7-iron in copper are 



found to be antiferromagnetic with the Neel temperature 
between 46 and 67 K 11 13 . 

The large Wilson-Sommefeld ratio Rw = 
(7r 2 /c|x)/(3/iB7) ~ 25, which can be deduced from 
the magnetic 10 and calorimetric measurements 14 , points 
to the presence of strong ferromagnetic fluctuations 
and indicates that the (antiferro) magnetism in 7-iron 
is likely to be frustrated by the competing magnetic 
fluctuations. Recent analysis 15 within the ab-initio spin 
spiral approach have also shown presence of long-range 
competing exchange interactions. Such an analysis 
assumes however the existence of the local magnetic 
moments and at the same time neglects electronic 
correlations. 

In the present paper we perform LDA and 
LDA+DMFT calculations to investigate the origin of 
weak antiferromagnetism of 7-iron, dominating types of 
magnetic fluctuations and possibility of the local moment 
formation in this substance. 



II. SPECTRAL PROPERTIES 

We first consider the results obtained for 7-iron in LDA 
approximation. 7-iron crystallizes in a stable face cen- 
tered cubic structure in temperature interval from 1183 K 
to 1667 K and it has the lattice parameter a = 3.656 A at 
1183 K 2 ' 3 . Band structure calculations have been carried 
out in LDA approximation 16 within tight-binding lin- 
ear muffin-tin orbital atomic spheres approximation (TB- 
LMTO-ASA) framework 17 . The von Barth-Hedin local 
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FIG. 1: (Color online) Left panel: the fatbands for ti g and 
e g orbit als in light (green) and dark (red) gray colors, re- 
spectively. Fatness corresponds to appropriate partial orbital 
contribution. Right panel: Iron density of states. Total DOS 
is shown by solid (black) line. Partial DOSes for t2 g , e g and 
sum of s + p orbitals are shown by light (green) , dark (red) 
and dashed-dark (blue) gray lines, respectively. 



exchange-correlation potential has been used 18 . Primi- 
tive reciprocal translation vectors have been discretized 
into 12 points along each direction which leads to 72 k- 
points in irreducible part of the Brillouin zone. 

The band structure together with the density of states 
are presented in the Fig. 1. On the left part of the figure 
the fatbands for the ti g and e g orbitals are shown by 
green and red colors, respectively (light and dark gray 
in the black-and-white version). The fatness coincides 
with the contribution of the corresponding partial DOSes 
shown on the right part of the Fig. 1. The bands of 
t2 g and e g symmetries hybridize in the vicinity of the L 
point and in K-T direction. In other symmetry directions 
the t2 g and e g manifolds hybridize weakly with s and p 
bands which span energy range from -8 eV to far above 
Fermi level (corresponding to zero energy). The t2 g states 
have a very flat region along X-W-L-K directions that 
is reflected in the DOS peak at 0.7 eV. At the Fermi level 
the partial t^g DOS has a deep. Another large peaks of 
the t2 g DOS are located at —1.3 and —2.6 eV. 

Although the e g partial DOS has a bandwidth almost 
equal to the t2 g counterpart, its shape is very different. 
The corresponding dispersion has a flat part at small neg- 
ative energy near T point (extended van Hove singularity, 
cf. Ref. 19), which results in the large peak of DOS just 
below the Fermi level at about —0.2 eV. This is in con- 
trast to a- iron 8 ' 20 , where peak of e g density of states is 
located very close to the Fermi level. As it will be shown 
below, this shift is of crucial importance for the mag- 
netic properties difference between a- and 7-iron. The 
smaller peak of the corresponding partial DOS is located 
at -3.4 eV. 

The Fermi surface obtained within LDA is shown in the 
Fig. 2. The four sheets that satisfy the equation for the 




FIG. 2: (Color online) 7-iron Fermi surface sheets. The 
colorcoding reflects contribution of the orbital states. The 
(red, green, blue) scheme is used for the color definition of the 
point where red is for e g , green is for t2 g and blue is for s + p 
orbitals, respectively. 



Fermi surface, £\^ F = 0, are colored such that amount of 
the appropriate color corresponds to the weight of par- 
tial contribution (we use the same colors as in Fig. 1: 
red for e g states, green for t2 g states and blue for s + p 
orbitals, respectively). The sheet a of the Fermi sur- 
face (Fig. 2a) is of mostly s, p, and e g orbital charac- 
ters. The sheets b and c (Figs. 2b, c) are mixture of t<i g 
and e g characters. The last sheet d (Fig. 2d) consists 
mostly of e g states. One should note that b and c sheets 
touch each other at the wavevector (0.57,0,0)27r and thus 
lead to the three bands crossing the Fermi level along 
T-X direction (see Fig. 1). Near the touch point these 
sheets have a cross-like features perpendicular to [0,0,1] 
direction with the small opposite incurvature produced 
by mostly t2 g states. This results in the approximate in- 
terband nesting of these crossed parts with close to zero 
wavevector and the intraband nesting with the wavevec- 
tor qA = (0.86, 0, 0)27r. The d sheet reminds the cube 
stretched along diagonals and it has also the cross-like 
feature. Its existence allows one to consider two addi- 
tional candidates for nesting vectors: within this sheet 
with qe = (0.48, 0, 0)27r and the vector connecting the 
sheets b,c and d, qc = (0.81, 0, 0)27r. 

In order to take into account correlation effects in 
3d shell of 7-iron we apply the LDA+DMFT method 
(for detailed description of the computation scheme see 
Refs. 21,22). The Coulomb interaction parameter value, 
U = 2.3 eV, and the Hund's parameter, I = 0.9 eV, used 
in our work are the same as in earlier LDA+DMFT cal- 
culations by Lichtenstein et al 23 for a-iron. The effec- 
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FIG. 3: (Color online) The imaginary parts of self-energies 
for t2 g (green) and e g states (red) plotted on the Matsubara 
energy grid for different temperatures (T=1290 K — circles 
and T=1550 K — squares). 



tive impurity model for DMFT was solved by quantum- 
Monte-Carlo (QMC) method with the Hirsch-Fye algo- 
rithm 24 . Calculations were performed for the value of 
temperature T « 1290 K which is just above the a- 
7 structural transition temperature. Inverse tempera- 
ture interval < r < f3 = l/fc#T was divided in 100 
slices. Four million QMC measurements were used in 
self-consistency loop within LDA+DMFT scheme and up 
to twelve million to refine data for spectral functions cal- 
culation with maximum entropy method 25 . 

The imaginary parts of self-energies for two values of 
temperature are presented in the Fig. 3. At low energies 
the behavior of the ^£(zco> n ) is qualitatively similar for 
the tig and e g orbitals. One can clearly see that increase 
of temperature does not change the frequency depen- 
dence qualitatively. The effective mass stays close to the 
bare value, m*/ra=1.2, and decreases slightly in temper- 
ature interval 1220 K < T < 1550 K, where 7- iron exists 
in nature. One should note that the obtained imaginary 
part of e g self-energy in 7-Fe has a quasiparticle-like 
frequency dependence, in stark contrast with the non- 
quasiparticle frequency dependence in a-phase 8 . This 
difference could be related to the peculiarities of the bare 
partial density of states near the Fermi level: the peak 
of partial DOS of e g states in 7-iron is shifted from the 
Fermi level in comparison with a- iron down to -0.2 eV 
and the states at the Fermi energy lie at the slope of peak. 
The very similar situation is observed in iron-based su- 
perconductors 26 . 

Partial densities of states obtained in paramagnetic 
LDA+DMFT calculation for ti g and e g electrons are 
presented in Fig. 4. The LDA+DMFT densities of states 
are slightly narrower than the LDA counterparts imply- 
ing weak correlation effects. This is in agreement with 
the small mass renormalization. One can observe that 
peak of e g density of states obtained in LDA approach 
is broadened in LDA+DMFT calculation. The shape of 
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FIG. 4: (Color online) The iron ti g (top panel) and e g 
(bottom panel) partial density of states obtained within 
LDA+DMFT method (solid lines) compared with bare LDA 
PDOS (filled area). 



tig density of states in LDA+DMFT approach resembles 
the LDA result with smearing of the peaky structures. 



III. MAGNETIC PROPERTIES 

To gain insight into favorability of different types of 
magnetic order in 7-iron, we analyze the momentum, q, 
dependence of generalized static magnetic susceptibility 
X q within LDA and LDA+DMFT approximations. The 
static magnetic susceptibility without correlation effects 
can be obtained as 
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where the Green function 



^ DA (^n) 



(io; n - H k 



, [i is the chemical potential and is the LDA- 
constructed Hamiltonian. Note that the temperature in 
Eq. (1) is introduced via the Fermi distribution function 
only. To analyze the contribution of different orbitals to 
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FIG. 5: (Color online) Contributions of different orbitals 
to magnetic susceptibility calculated according to the Eq. 
(5) in high symmetry directions. Black line corresponds to 



Xq d - Red, green and blue lines show x<t 



o,t 2 

Xq 



and Xq 6s 2s ? respectively. Solid (dashed) lines correspond 
to T=0 (T=1290 K). Wavevectors are measured in units of 
2tt. 



the susceptibility, we represent Green function 
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where {\m)} is an orbital (LMTO) basis and ip£ m (e a k) 
are LDA eigenvectors (eigenvalues) written in orbital rep- 
resentation (a is a band index). In this notation the 
equation (1) can be rewritten as 
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For the following analysis we split the susceptibility as 

(4) 
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where the first term, x q ' d , corresponds to restricting the 
mi ; 2 sum in Eq. (3) over d-orbitals only, while the second 
term, Xq' rest ? contains the rest. This separation can be 
repeated further to analyze contributions to Xq d from 
different orbitals: 
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The results of such decomposition are presented in the 
Fig. 5. The intra-orbital contributions to the suscepti- 
bility at zero temperature, Xq 6s &E an d Xq t2s t2s ? are °f 
the same magnitude and varying in "counter-phase" and 
thus compensating to large extent the q dependence of 
each other. The e 9 — e g contribution has weak momen- 
tum dependence, favoring ferromagnetic ordering (maxi- 
mum at T point) and features at the nesting wavevectors 



FIG. 6: (Color online) Contributions of different orbitals to 
irreducible susceptibility calculated according to the Eq. (1 
(dashed lines) and in LDA+DMFT approach (Eq. (7), solid 
lines) in high symmetry directions at T =1290 K. Colorcoding 
and units repeats the previous picture. 



qe and qc, discussed in Sec. II. Note that the mag- 
netic behavior of e g orbitals is very sensitive to the tem- 
perature, which is related to the peculiarities of the e g 
band dispersion (the flatness of the spectrum along the 
direction T— L and the small size of the d sheet of the 
Fermi surface, see Fig. 2d). The momentum dependence 
of t2 g — t2 g contribution is almost negligible and has a 
maximum at the wavevector qx = (27r,0,0), related to 
the intraband nesting of the c Fermi surface sheet. The 
largest part of the momentum dependence of suscepti- 
bility comes from e g — t^g contribution, which, at zero 
temperature, is maximal at the point qA and negative 
and large by magnitude in the vicinity of q = point 
due to small momentum transfer between electron-like 
(mainly t2 g -derived) Fermi-surface sheet c and hole-like 
(mainly e^-derived) sheet b. The maximum of the result- 
ing susceptibility x q ' d is obtained close to the wavevector 
qA- Increasing temperature shifts maximum of e g — t2 g 
contribution, together with the maximum of the <i-orbital 
susceptibility to the wave vector qx, stabilizing even fur- 
ther the antiferromagnetic fluctuations. The wave vector 
qx of the maximum of susceptibility corresponds to an- 
tiferromagnetic structure with alternating orientation of 
magnetic moments in adjacent layers of face-centered cu- 
bic crystal structure. This magnetic configuration was in 
fact observed in the precipitates of 7- iron 27 . 

The effects of electron-electron interaction can be 
treated within LDA+DMFT approach. Since, in gen- 
eral, interaction produces vertex correction to a single 
bubble considered above, we neglect for sake of simplic- 
ity the frequency dependence of these vertex corrections, 
introducing the frequency-independent vertex r irr , such 
that 
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H(icj n ) is DMFT self-energy with subtracted double 
counting term, Vd is a projector onto <i-orbitals and 5ji is 
a change of the chemical potential in DMFT with respect 
to LDA value. 

Orbitally-resolved contributions to the irreducible sus- 
ceptibility in LDA+DMFT approach are presented in the 
Fig. 6. One can see that the DMFT self-energy correc- 
tions lead to suppression of irreducible susceptibility, the 
momentum dependence becomes more flat in compari- 
son to LDA, except the dip in the vicinity of the point 
q = 0, which appear due to small momentum transfer 
processes between sheets b and c of the Fermi surface. 
The flat region implies close competition of the antifer- 
romagnetic fluctuations with the wavevectors along the 
directions qx — (7r,7r,7r) — (37r/2, 37r/2, 0). According to 
the general ideas of spin- fluctuation theory 28 , the weak 
momentum dependence of the irreducible susceptibilty 
can be also attributed to the partial presence of local 
moments, discussed below. 

To get further insight into the interplay of different 
magnetic fluctuations in 7-iron, we consider the uniform 
magnetic susceptibility; the latter can give a key for 
understanding the role of magnetic fluctuations. The 
uniform magnetic susceptibility x(T) in the paramag- 
netic state of 7-iron was extracted from the LDA+DMFT 
simulations as a ratio of the induced magnetic moment 
by a small external magnetic field and the field mag- 
nitude 26 ' 29 . The temperature dependence of x _1 (^) is 
presented on Fig. 7. We note the absence of fulfillment 
of the Curie- Weiss law 



X(T) 



Mew 



3(r-0cw)' 
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up to highest considered temperatures and the minimum 
of x _1 ( T ) at T * - K - The effective magnetic mo- 
ment, extracted from the slope of the inverse susceptibil- 
ity in the low-temperature region (1200-1550 K), /icw 
= 5.75/iB, is close to the experimentally observed value, 
Mew = 7.47//B 9 ' 10 - On the other hand, despite the Curie- 
Weiss law is not satisfied, roughly estimating the effective 
magnetic moment value, /icw, from high- temperature re- 
gion (2500-4000 K) we find smaller value /icw ~ 4/iB 
(which approximately corresponds to spin S = 3/2). 

In order to analyze the role of peculiarities of band 
structure on non-monotonous temperature behavior of 
x(T), we calculate XqLo(^) projected onto pair sets of 
orbitals as in Eq. 5. The results are shown in Fig. 8 
and inset of the Fig. 7. The overall temperature depen- 
dence of XqLo(^) repeats that of x(T), being however 
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7: (Color online) Temperature dependence of the 
inverse uniform magnetic susceptibility calculated within 
LDA+DMFT (blue circles) and experimental data 9 (red cir- 
cles), red line corresponds to the least square fit to Curie- 
Weiss law. Shadow covers the temperature range of 7-phase 
existence. Inset shows the inverse total (black) and orbital 
(red-e^, green-fog) contributions to Xq=o- 
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FIG. 8: (Color online) Temperature dependence of Xq=o cal- 
culated within LDA+DMFT. Top panel - Xq=o> middle panel 
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substantially weaker. The fo^ contribution to Xqlo nas 
maximum at T ~ 2000 K, which originates from the max- 
imum of ^-projected DOS above the Fermi level at the 
energy about 0.3 eV, similarly to pnictides 26 . The e g - 
t<i g contribution behaves similarly by magnitude at low 
enough temperature (T < 1000 K) but has a negative 
sign. The contribution of e g orbitals decreases almost 
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linearly with increasing temperature. This is connected 
with strong (relatively to t^g orbitals) correlated char- 
acter of e g orbitals. Such distinct behavior of different 
orbitals contributions results in the shift of maximum of 
total (i-orbital irreducible susceptibility to approximately 
the temperature T*, making it close to the position of 
uniform susceptibility maximum. We therefore interpret 
T* as a crossover temperature which corresponds to the 
change of the dominating magnetic fluctuations. The 
ratio of total uniform susceptibility and irreducible one 
(Stoner enhancement factor) at T ~ 1290 K is about 10, 
which means that ferromagnetic fluctuations are strong 
in the temperature interval in the vicinity of T*. Such a 
large ratio also explains strong temperature dependence 
of x(T) in comparison with x q lo(^)- 

To investigate the possibility of the local moment for- 
mation in 7-iron the real part of the dynamic local mag- 
netic susceptibilities for different temperatures have been 
calculated. In Fig. 9 we present these functions for 
different temperatures, rescaling both the susceptibility 
and frequency by temperature. For comparison, we also 
present on the inset the corresponding result for a-iron 
(see also Ref. 8). The dynamical mean- field theory, which 
neglects intersite magnetic exchange and therefore has no 
other low-energy scales apart from frequency and temper- 
ature, is expected to yield the local magnetic susceptibil- 
ity, Xioc(^) = {^-/T)f{oo/T), for the system with the local 
moments (where f(x) — >• at x — >• oo). Such a depen- 
dence also naturally provides the static nature of a single 
spin, xioc 5(uj) at T — » 0. One can see that the suscep- 
tibility for a-iron obeys fully the above discussed scaling 
form, while for 7-iron some deviations from it are ob- 
served. The peak widths of the 7-iron susceptibilities are 
larger than those for a-iron, and the height of the peak 
(in units of I/Q^bT)) varies with temperature. There- 
fore the local moments in 7-iron are not well defined, 
which is also confirmed by the quasiparticle form of the 
self-energy obtained in Sect. II. However, the shape of 
the dependence Tx(co>/T) changes only weakly with tem- 
perature allowing to estimate the size of the "effective" 
local moments from the values of static local suscepti- 
bility, /i 1 2 oc = 3Txioc(0). We obtain ji\ oc « 3/iB in the 
temperature region 1200K-1400 K. 

To estimate exchange interactions we perform the 
mapping of the considered electronic system to the ef- 
fective Heisenberg model. Due to presence of different 
competing magnetic orders we consider a rough way to 
extract the exchange integrals using the electronic prop- 
erties in the paramagnetic phase at finite temperature. 
To this end we compare a mean-field-like momentum de- 
pendence of the static magnetic susceptibility, Xq? within 
the effective Heisenberg model with exchange parameters 
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(where x° = X°CO = mL/( 3T ))> with the E <1- ( 6 )> which 




FIG. 9: (Color online) Local magnetic susceptibility of 7-iron 
for different temperatures. The inset shows the results for 
a-iron. 



yields 



Jq = -4/4 (Xq r ) X + const. 



(11) 



Using the results for Xq r within the LDA+DMFT 
method one can obtain the constant in the Eq. (11) if one 



fix that J q = 0. At T 

-2380 K and max q J q 



1290 K we obtain min q J q 



q=0 : 

while for T 
max q J q 



Ja 



Ja 



580 K min q J q 



Ja 



q=qx = 1172 K, 
q =o = -5549 K and 
q = qx — 1442 K. Therefore, the exchange in- 
tegrals with wavevectors near qx vary slightly with tem- 
perature, while the exchange J q =o acquires strong tem- 
perature dependence, which implies that the local mo- 
ment picture is better defined for magnetic states with 
the wavevectors near qx- 

To estimate roughly the Neel temperature we consider 
the solution of mean-field-type equation 



T MF _ 
^Neel — 



ioc^q= 



■<ix)t=t£ 



3(2^ B ) 2 



(12) 



performing linear interpolation for J q = qx (^) an d ^\ oc (T) 
dependence. From this equation we obtain = 793 

K. The Neel temperature is therefore strongly overesti- 
mated within our approach with respect to the experi- 
mental value. This agrees with the presence of strong 
frustration in 7-iron; the result (12) should be consid- 
ered as an estimate for characteric temperature of enter- 
ing antiferromagnetic fluctuations-driven regime, which 
is of the same order of magnitude as T* obtained above. 
Although the competition of different magnetic orders 
takes place already in the Heisenberg model on the face- 
centered cubic lattice 30 (making the mean-field treat- 
ment inadequate), it is expected to yield even stronger 
reduction of TNeel in itinerant system due to strong sup- 
pression of sublattice magnetization by paramagnons, 
having small gaps in their spectrum in the symmetric 
Brilloun zone points, different from q = qx- 
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Considering analogously to weak itinerant ferromag- 
nets the Rhodes- Wohlfarth ratio 28 , R = Pcw/Ps, where 
Pew ~ 3.13 is effective moment extracted from the Curie- 
Weiss law and p s = 0.75 is a saturation moment, we 
obtain a large value R ~ 4.17, which is typical for itiner- 
ant but not the local moment magnets. This shows that 
in general 7-iron can be classified as a weak itinerant 
antiferromagnet with strong influence of the frustration 
effects. 



IV. CONCLUSIONS 

We have considered the electronic and magnetic prop- 
erties of paramagnetic 7-iron. The shift of the DOS peak 
below the Fermi level in 7-iron causes the dramatic differ- 
ence in the electronic and magnetic properties between a- 
and 7-iron. The position of this peak is therefore crucial 
for understanding the magnetic properties which is sim- 
ilar to recent study of pnictides 26 . Other band structure 
features appear to be important for the antiferromag- 
netism of 7-iron, in particular, the e g -t2 g hybridization 
favors to the nesting of Fermi surface. The obtained anti- 
ferromagnetic state with the wavevector close to (27r, 0, 0) 
is found to compete strongly with the other incommensu- 
rate spin-density wave instabilities. Observed tendency 



to the magnetic frustration can explain the small Neel 
temperature of 7-iron. 

The account of correlation effects in 7-iron allows one 
to conclude that: i) the correlation strength is from mod- 
erate to weak; ii) the local moments are not well formed 
albeit the "effective" local moments can be introduced 
with /ii oc & 3/iB in the temperature range 1200K-1400 K. 
In general, 7-iron should be considered as a weak itiner- 
ant antiferromagnet with some importance of correlation 
effects. 
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